Recovering the water-wave profile from pressure measurements 
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P_( Abstract 

^ A new method is proposed to recover the water-wave surface elevation from pressure data 

obtained at the bottom of the fluid. The new method requires the numerical solution of a 

1 " I nonlocal nonlinear equation relating the pressure and the surface elevation which is obtained 

from the Euler formulation of the water-wave problem without approximation. From this new 
equation, a variety of different asymptotic formulas are derived. The nonlocal equation and 
the asymptotic formulas are compared with both numerical data and physical experiments. 
The solvability properties of the nonlocal equation are rigorously analyzed using the Implicit 
Function Theorem. 

^ 1 Introduction 

In field experiments, the elevation of a surface water-wave in shallow water is often determined by 
J> measuring the pressure along the bottom of the fluid, see e.g. [4], [6], [l8], [H], [21], [25]. A variety 

^ of approaches are used for this. The two most commonly used are the hydrostatic approximation 

and the transfer function approach. For the hydrostatic approximation [9l I17j. 

, , P(x,—h,t) , , , 

P9 

where g is the acceleration due to gravity, h represents the average depth of the fluid, p is the fluid 
density, P{x, —h, t) is the pressure as a function of space x and time t evaluated at the bottom of 
the fluid z = —h, and r]{x,t) is the zero-average surface elevation. Throughout, we assume that 
all wave motion is one-dimensional with only one horizontal spatial variable x. The hydrostatic 
approximation is used, for instance, in open-ocean buoys employed for tsunami detection, see |21j . 
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The transfer function approach uses a hnear relationship between the Fourier transforms T of 
the dynamical part of the pressure and the elevation of the surface [9l HSj [HI [IT] : 

T{7](x,i)) (k) = coshikh)T{pix,t)/g} (k), (2) 

where p{x, t) = {P{x, —h, t) — pgh) / p is the dynamic (or non-static) part of the pressure P{x, z, t) 
evaluated at the bottom of the fluid z = —h, scaled by the fluid density p. In this relationship, rj 
and p are regarded as functions of the spatial coordinate x, with parametric dependence on time t. 
It is equally useful to let t vary for fixed x, as would be appropriate for a time series measurement, 
which results in extra factors of the wave speed c{k), due to the presence of a temporal instead of 
a spatial Fourier transform. 

It is well understood that nonlinear effects play a significant role when reconstructing the surface 
elevation for shallow- water waves or for waves in the surf zone (see [5l |6l [25], for instance). Since 
nonlinear effects are not captured by the linear transfer function ([2]), different modifications of 
([2]) have been proposed. One approach is to modify the transfer function to incorporate extra 
parameters {e.g., multiplicative factors, width scalings) that are tuned to fit data [I8l [HI [25] . A 
less empirical approach is followed in [5] and [20] where corrections to the transfer function are 
proposed based on higher-order Stokes expansions. Bishop & Donelan [S] examine the empirical 
approaches and argue that the inclusion of the proposed parameters is not necessary as errors from 
inaccuracy in instrumentation and analysis are likely to outweigh the benefit of their presence. We 
do not include any of the modified transfer function approaches in our comparisons. 

Bishop and Donelan [6] acknowledge that the linear response cannot accurately capture nonlin- 
ear effects. While both the hydrostatic model and the transfer function approach are accurate on 
some scales, they fail to reconstruct the surface elevation accurately in the case of large- amplitude 
waves, as might be expected. Errors of 15% or more are common, as is shown and discussed below. 
In order to address the inaccuracies of the linear models, nonlinear methods are required. With the 
exception of recent work by Escher &: Schlurmann ^_3j and Constantin &: Strauss ^ , few nonlinear 
results are found in the literature. Escher &: Schlurmann [13] provide a consistent derivation of ([2]) 
and offer some thoughts about the impact of nonlinear effects. Starting from a traveling wave as- 
sumption, Constantin & Strauss [8] obtain different properties and bounds relating the pressure and 
surface elevation. However, they do not present a reconstruction method to accurately determine 
one function in terms of the other. 

One way to obtain an improved pressure-to-surface elevation map is to use perturbation methods 
to determine nonlinear correction terms to ([2]). Several such approaches are given below, and we 
include them when comparing the different methods. Our main focus, however, is the presentation 
of a new nonlocal nonlinear relationship between the pressure at the bottom of the fluid, and the 
elevation of a traveling-wave surface that captures the full nonlinearity of Euler's Equations. The 
advantage of this approach is that 

1. it allows for the surface to be reconstructed numerically from any given pressure data for a 
traveling wave, 

2. it provides an environment for the direct analysis of the relationship between all physically 
relevant parameters such as depth and wave speed, 

3. and it allows for the quick derivation of perturbation expansions such as the ones mentioned 
above. 
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Figure 1: The fluid domain D for the water wave problem. An ideahzed pressure sensor is indicated 
at the bottom. In our calculations the pressure measurement is assumed to be a point measurement. 



4. Although our approach is formally limited to traveling waves, it can be applied with great 
success to more general wave profiles that are not merely traveling. This is illustrated and 
discussed below. 

In what follows, we derive these nonlocal relations and demonstrate their practicality. We 
compare results from the nonlocal formulation with those from the linear approaches and different 
nonlinear perturbative models, using both numerical data for traveling waves in shallow water, and 
experimental data obtained at Penn State's Pritchard Fluid Mechanics Laboratory. We demon- 
strate the superiority of the nonlocal reconstruction formula for a large range of amplitudes. In 
addition, using the Implicit Function Theorem, we analyze the nonlocal formulation in order to 
demonstrate its solvability for the surface elevation given the pressure. 

2 A nonlocal formula relating pressure and surface elevation 

Consider Euler's equations describing the dynamics of the surface of an ideal fluid in two dimensions 
(with a one-dimensional surface): 



(t>xx + ^zz = 0, (x, z) G D, (3) 

^z = 0, z = -h, (4) 

ilt + r]x<j)x = 4>z, z = r]{x,t), (5) 

(l>t + l{(t}l + 4) +9ri = 0, z = r,{x,t), (6) 



where (p{x, z, t) represents the velocity potential of the fluid with surface elevation r/(x, t). As posed, 
the equations require the solution of Laplace's equation inside the fluid domain D, see Figure [ij 
If the problem is posed on the whole line x G M, we require that all quantities approach zero at 
infinity. If periodic boundary conditions are used then all quantities at the right end of the fiuid 
domain are equal to those at the left end. Below we work with the whole line problem, stating only 
the results for the periodic case. 

Following [26j, let q{x,t) represent the velocity potential at the surface z = 7]{x,t), so that 

q{x,t) = (j){x,r](x,t),t). (7) 
Combining the above with equation ([5|, we have 
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(t>z = 'rit + {Qx - 4'z'nx)r]x, 

obtained by taking an x-derivative of ([T]). This allows us to solve directly for (f>z in terms of i] and 

q- 

, Vt+ VxQx 

Using ([5| again gives an expression for (j)^, while taking a t-derivative of ([T]) leads to an expression 
for (pf- 

^-=^J-^, ^, = ,,- ^^^]^^Y'\ (9) 



Substituting these expressions into the dynamic boundary condition ^ we find 



Qt + i^1x+9r]-i^ = 0, (10) 



1„2 , _ \ (j]t+JlxVxf 
2 l + r/2. 

after some simplification. 

Next, we restrict to the case of a traveling wave moving with velocity c. We introduce = x — ct, 
so that X and t-derivatives become ^-derivatives, the latter multiplied by — c. The Bernoulli equation 
(10) becomes a quadratic equation in q^: 

2 ^ 1 "lli^i 

U +gV- 7^ -^^rr- ^ 



1 o 1 viiQi -cY 
c., + ^^| + ..-^^Y^ = o. (11) 



Solving this quadratic equation, we find 

g5 = c±^(c2-2OT)(l + 7?|), (12) 

where for c > we choose the — sign, to ensure that the local horizontal velocity is less than the 
wave speed [H]. Similarly, for c < the + sign should be chosen. 
Substituting this result into ([s]) and Q, we find 



where we have chosen c > 0, without loss of generality. This simple calculations allows us to to 
express the gradient of the velocity potential at the surface directly in terms of the surface elevation. 

Returning to the original coordinate system {x,z,t), let Q{x,t) = (p{x, —h,t), the velocity 
potential at the bottom of the fluid. Inside the fluid, we know that the Bernoulli equation holds: 

^t + l {€ + +9Z + ^^^^ = 0, -h<z< v{x, t). (14) 
z p 

Evaluating this equation at z = —h, we find 

Qt + lQl-gh+^^^^^ = o. (15) 

z p 
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Moving to a traveling coordinate frame as before, we obtain a quadratic equation for Q^. Solving 
for we find 



= c - - 2p, (16) 

where p(^) represents the non-static part of the pressure at the bottom in the traveling coordinate 
frame, scaled by the fluid density: p{^) = P{x — ct, —h, t)/p — gh. For consistency with our previous 
choice, we work with the — sign again. Next, we connect the information at the surface with that 
at the bottom of the fluid. 

Within the bulk of the fluid D, 

% + (pzz = 0, (17) 



where the boundary conditions given in (13) and (16) must also be satisfied. We can write the 
solution of this equation as 

1 f'°° 

</,(^, z) = — e**-'%(/fc) cosh {k {z + h)) dk, (18) 

where the boundary condition for 0^ at z = —h is satisfied. For the boundary condition at the 
bottom for (ht we find 



27r 

so that 



1 

- / ike'''^^{k) dk = c- - 2p, (19) 



iA;^(A:) = 2^C(5(fc) - J^|x/c^^2p} (A;), (20) 

where 5{k) is the Dirac delta function and denotes the Fourier transform: J^{y{(^)}{k) = 
I^oo y(^^ exp(— iA;^)d(^. Evaluating z) at the surface z = t], we have 

1 f°° 

(Pei^, V) = — / e'^^ik-^{k) cosh {k [r] + h)) dk 

= c- ^ J e^^^ cosh (A; (r/ + /i)) I yjc^ - 2p} (A;) dk. 



Using the boundary conditions given in (13), we find the nonlocal relationship 



— 2grj 1 
1 + r/| ~ 2^ 



/ e'^^ cosh {k{ri + h))T[ yjc^ - 2p\ {k) dk. (21) 



Equation (21) is the main result of this paper. It provides an implicit relationship between the 
surface elevation of a localized traveling wave ry(x) and the pressure measured at the bottom of the 
fluid p{^). In the rest of this paper we investigate how this relationship may be used to compute 
viO if P{0 is known, and how different asymptotic formulas may be derived from it. 
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Remarks. 



• In order to extend the above to periodic boundary conditions, we use the periodic generaUza- 
tion of the formulation of Ablowitz, Fokas, & Musshmani AFM, see [llj . FoUowing the steps 
outHned above, this aUows for the derivation of a relation between the surface elevation of a 
periodic traveling wave and the pressure at the bottom: 



2gr] 



1+7?| 



1 

2^ 



J2 e'''^ cosh {k{rj + h))Pk 



(22) 



where Pk 



r27r _ 

Jo e 



2p{^) d^. In what follows, we will use either (21) or (22). 



In the above, we have assumed there exists a smooth solution to the water wave problem 
(ppl). Given a speed c and a non- hydrostatic pressure profile p as inputs, we aim to solve (21 ) 



for rj. However, these inputs cannot be independent of each other: indeed arbitrary pressure 
profiles will not lead to surface elevations corresponding to solutions of ([3]j6]). One expects 
that for a given speed c, there exists a unique surface elevation rj and associated pressure 
profile p. In order to back up this intuition, we require another relation between rj, p and 
c. Such a relation is found by taking a derivative with respect to 2 of (18) and equating the 



result to the right-hand side of the second equality in (13). Finally, (20) is used, resulting in 



2gr] 



—I 
2^ 



Akx 



sinh(A;(r? + /i))J"|c- ^/c^ - 2p^{k)dk. (23) 



The system of equations (21) and (23) may be solved to obtain both 77 and p, given c. We 



will not pursue this issue further and content ourselves with establishing a map from p to r]. 

For the purposes of the question considered in this paper, the above is not an issue: we assume 
that the given pressure originates from experimental observations and hence corresponds to 
the unique solution of (3][6), to the extent that the Euler equations provide an accurate model 
for the water wave problem. 

Obtaining the pressure at the bottom from the surface elevation. An explicit 
nonlinear relationship for the pressure at the bottom in terms of the surface elevation may 
be obtained directly from the approach of Ablowitz, Fokas and Musslimani. Consider the 
relationship (the one-dimensional version of Equation 1.11 in [Ij) 



J — ( 



^ikx±\k\{ri+h) 



e ' '~' {sgn{k)rit - i Qx) dx = -i 



X, —h, t)dx, k G 



^0- 



(24) 



Changing to a moving frame of reference and substituting the expressions ( |12|16D into this 
relation one obtains 



Jkx±\k\{rj+h) 



icsgn{k)rix + ( c - y/ (c^ - 2gr]){l + r]"^] 



dx 



ikx 



(c- \/c 



2pj dx. 
(25) 
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The right-hand side is essentially the Fourier transform of the quantity c— sj <? — 2p. Inverting 
this transform, one may solve for the pressure at the bottom in terms of the surface elevation. 
Equation (25) provides an alternative to (21) and can be used in its stead. The formula 



(25) is advantageous if one wishes to compute the pressure, given the surface elevation. The 



presence of a small-divisor problem (seen by linearizing the left-hand side integrand about 
77 = 0) when c is near its shallow- water limit value ^/gh indicates that (21) is to be preferred 
over ( |25[ ) to reconstruct the elevation t/, given the bottom pressure p. 

In [8] , Constantin and Strauss derive various properties of the pressure underneath a traveling 
wave. It should be possible to re-derive these properties directly from (21) or (22). This is 
not pursued further here. 



3 Existence and uniqueness of solutions to the nonlinear formula 



In this section, we analyze (21). Among other results, using the Implicit Function Theorem, we 



show that the nonlocal relation (21 ) gives rise to a well-defined map from the pressure profile to the 



surface elevation: given the pressure profile p at the bottom, (21 ) defines a unique surface elevation 
r/. In other words, we can expect the asymptotic and numerical methods employed in the next 
sections to produce faithful approximations to the true solution. 
Define the operator F, parameterized by c E M, by 



F{r],p) 



2grj 



1 + ^^ 



1 

2^ 



/oo 
e'^'"" cosh(A:(r/ + h))F\c- - 2p \ {k)dk. 
-00 ^ ^ 



(26) 



Note that F{ri,p) = is equivalent to (21). Using the Implicit Function Theorem, we wish to show 
that the equation F{r],p) = has a solution profile rj, given sufficiently small pressure p. We have 
that F{0, 0) = 0. In order to apply the Implicit Function Theorem, we need to define appropriate 
Banach spaces for which the operator F is defined. First we seek a suitable space for rj. An obvious 
choice is r/ G C^[M,]R], i.e. ry is a continuously differentiable function which vanishes at infinity. 
This space is supplied with the usual norm: 



\v\\c^ = sup \r]{x)\ + sup |r/'(x)|. 



(27) 



If ll^/llci < (? /2g then c — (c^ — 2gri) /(I + r/^) represents a continuous function of x. Hence 
we are motivated to define the image of F in C[M,M]. Consequently, we wish for 



e'^^ cosh{k{r] + h))F |c - ^/c'^ - 2pYk)dk, 



to be a continuous function of x. For finite H^/Uci, this nonlocal term is a continuous function of x 
if 



/oo 
cosh{k{Mci + h)) f\c- Vc'^ -'^p\{k) 



dk < 00, 



(28) 



and if the integrand of the nonlocal term is a continuous function of x for every k (see Theorem 2.27 
on pg. 56 of [15j). Let us consider the second condition, namely the continuity of the integrand. 
Since by assumption rj is continuous, the continuity in x of the integrand requires 
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An application of the Cauchy-Schwarz inequality gives 



sup 

k 



< oo. 



(29) 



{c-Vc2-2p}(fe) 



< 



1 



-oo 1 ~l~ 1-^1 

Hence, we impose the following condition on the pressure p: 

2 



(1 + C - \/c2 - 2p 



dx 



/oo 
(1 + c- \/c2 - 2p 
-oo 



dx < oo. 



(30) 



Next, we return to the first condition (|28|). Due to the presence of the hyperbolic cosine, 



we expect that it is necessary for J^{c — sj (? — 2p}{k) to have sufficient decay for large |fc|. Let 
M > h + \\r]\\ci- Starting from the integral in ( p8| ), we apply the Cauchy-Schwarz inequality again 
to find 



r cosh(fe(||7?||ci +/i))e-*^l'=le^^l'=l T ic - ^/c^ - 2p\{k) 

J — oo ^ ^ 



dk 



< 



< C 



cos 



MHWvWc^ + h))e-^'\''\) dk 



,2M\k\ 



{c - - 2p}(fc) 



2 \ 2 

dk\ 



„2M|fc| 



v/c2 - 2p^(k) 



2 \ 2 

dk\ 



for some constant C . Thus, if the conditions 



r ^2M\k\ ^|c-y^^3^|(A:) 

J — oo 



and 



(1 + \xf) c- Vc2 - 2p 



dk < oo 



dx < oo 



hold, the function F{r],p){x) is continuous for 77 G C^. 



(31) 
(32) 



Having found the conditions (31 32), we now determine an appropriate function space for p so 



that they are satisfied. The following theorem due to Paley and Wiener (Theorem 4 pg. 7 of [23j) 
is helpful. 

Theorem 3.1 If w{z) (with z = x + iy) is an analytic function in the strip —A < y < where 
X, fi > and 



\w{x + iy)\'^dx < 00, —A < y < fj., 
then there exists a measurable function w{k) such that 
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\w{k)\^e^^'^dk < oo, / \w{k)\^e-'^^^dk < oo, 



and 

w{x + iy) = lim / —w{k)e'^'~''+''yUk, -X<y<n 

where the limit is to be understood in the mean-square sense. 

In other words, the Fourier transform of w{x) exists and it has decay as specified above. In 
particular, for any M < mm{X, fi} 

e^^'^''^\wik)fdk= / e^^"'\w{k)fdk+ / e-^^"'\w{k)\^dk, 



oo JO 



oo 

J — oo 

rO 

< I e^f'^\w{k)\^dk+ / e'^^''\wik)\'^dk <oo. 











The theorem implies that a sufficient condition for (31) to hold is that c — — 2p is an 
analytic function of z within a strip of width at least 2M centered around the real axis and that 
it is square-integrable along lines parallel to the real axis within this strip. Of course, the presence 
of the square root is a hinderance to the analyticity of the function. Consequently, we require 
that \p\ < everywhere in the strip. This condition also implies the square-integrability of the 
function if p is square-integrable. Indeed, the function 

f{z) = c - Vc^ - 2z, 

is analytic (and hence Lipschitz) in a neighborhood of the origin for which \z\ < (? j2. Hence, for 
all zi, Z2 in such a neighborhood of the origin we have 

|/(^1)-/(^2)|<C|Z1-Z2|, (33) 

for some constant C. In particular, since /(O) = 0, 

\f{z)\<C\z\, (34) 

uniformly for all l-z] < 5 < i.e., the constant C is independent of z. Next, consider the 

function f{p{z)), where p{z) is a function which is analytic and bounded in the strip of width 2M. 
This implies 

\f{p{z))\<C\p{z)\. (35) 

As above, the constant C is independent of p{z) and thus of z, provided \p{z)\ < Thus 

the square-integrability of p{z) implies the square-integrability of c — \J (? — 2p for |p| < 5 < 

for every z in the strip. Thus, if p(z) is an analytic function in the strip of width at least 2M, 
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square-integrable along lines parallel to the real axis, and bounded in the strip, the first condition 
(Mf holds. 

Another theorem due to Paley and Wiener (Theorem 2 pg 5 of [23]) allows us to bound \p\ in 
terms of the L^-norm. 

Theorem 3.2 Let w{z) be analytic in the strip —A <y<^ with /x, A > and 

\w {x + iy)\^dx < oo^ —X<y<fi, (36) 
then for any z in the interior of the region 



w{z) 



1 

271 



w{x + ifi) 



dx 



1 

2^ 



w{x — iX) 
X — iX — z 



dx. 



(37) 



In particular, an application of the Cauchy-Schwarz inequality shows that for any y G [— A + e, ^ — e] 
with e > 0, w(z) is bounded in terms of the norms of w{x + ip.) and w{x — iX). 

Collecting these ideas, we choose the pressure p to be in the space of analytic functions in the 
symmetric strip of width 2M about the real axis such that 



{I + \x\'^)\p{x + iy)\'^dx < oo, -M<y<M. 



(38) 



Note that this condition guarantees that the second condition (32) is also satisfied. Let Hm denote 
the space defined by (38). It is endowed with the norm 



WPWHm = sup 

\y\<M 



[1 + |xp)|p(x + iy)\^dx 



1/2 



(39) 



We claim that Hm is a Banach space. Indeed, with the obvious definitions of addition and scalar 
multiplication for elements p G Hm-, Hm is a vector space. It is straightforward to verify that (39) 
defines a norm. Thus, it remains to verify completeness. Let {fk\ be a Cauchy sequence in Hm- 
With the norm above, this sequence converges to a complex-valued function / defined on the strip 
since for fixed y, {fk] defines a Cauchy sequence in the space with weight (1 + |xp) and has the 
limit /(•, y) in this space for each y. We define the function 



/(x,y) = lim fk{x + iy), y fixed. 

fc— >oo 



(40) 



Since {fk} is a Cauchy sequence, for every e > there exists an such that for n,k > N 

\\fn-fk\\HM<^- (41) 



This implies 

/oo 
(1 + \x\'^)\fnix + iy) - fk{x + iy)\^dx < e, 
-oo 

for every \y\ < M. Letting — )• oo in the above integral we obtain 

(1 + \x\^)\fn{x + iy) - f{x,y)\^dx < e, 



(42) 



(43) 
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for every \y\ < M and thus fn^f^n the Hm norm. Theorem 3.2 imphes the pointwise bound 

\w{x + iy)\ < C||zi;(x + 



Ha 



(44) 



for x + iy in the interior of the strip. Consequently, convergence in Hm imphes uniform convergence 
on compact subsets of the strip. From Morera's Theorem, it foUows that the Cauchy sequence of 
analytic functions converges to an analytic function, thus the space Hm is complete. 

We now state the theorem for the existence of a map from the pressure beneath a traveling 
wave to the surface elevation of the wave. 

Theorem 3.3 Let p and rj be the bottom pressure and surface elevation, respectively, obtained by 
solving the Euler equations augmented with (14)- Assume that p £ Hm+e for some M > h, e > 
and that \\rj\\ci < min[M — h, (? /2g\. Then for fixed c ^ and sufficiently small p, the equation 



2gr] 



1 

2^ 



e'^"" cosh(fc(r/ + h))F |c - \/c2 - 2p}(fc)dA;, 



has a solution rj. Further, if p is the true pressure consistent with (3 6p^), then the only solution 
for rj is that which solves the stationary water wave problem with speed c. 



Proof. Let M > h and e > 0. By Theorem 3.2 there is a ball V around the origin in Hm+c, 
i.e. V = {p£ HM+e ■ \\p\\HM+e < '^}) such that 



sup \p{x + iy)\ < S < —. 

\y\<M+e/2 ^ 

By definition, p \s, bounded, analytic and square integrable along lines parallel to the real axis. 
Then the function c — \J (? — Ip is also analytic and square integrable along lines parallel to the 
real axis in a strip of width 2{M + e/2) symmetric with respect to the real axis. Using Theorem 



3.1 



J — oo 



dk < oo. 



For p £ y, following the discussion preceeding Theorem 3.1, define the reconstructed functions 

d)^'^ and (h^'^ as 



1 r°° ( ~i 

= — / (^^''co^\^(k{z^K))T\c- ^/^?^p\(k)dk, (45) 
cj)^''{x,z;p) = — J -ie'^'"' sinh(A;(z + h))T \c - ^ c^ - 2p| {k)dk. (46) 

Then ^^'^ and ^^'^ are harmonic (and thus smooth) for all x G M and |z + /i| < M. Let K = 
min (M - h,c^/2g) and define the ball U = {r] e : M\ci < R}. Then G : U x V C^[R,R], 
defined as 

G{ri, p) = -0^'^(x, r]-p) + ^ [(/-^'^ (x, 7?; p)Y + ^ [0^'"(x, r?; p)]^ + ^i], 

2c 2c c 
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is a continuously differentiable function with G{0, 0) = 0, as is readily verified by computing its 
second variation evaluated at {rj,p) = (0,0). The Frechet derivative of G with respect to rj at the 
origin is 

GJO,0)v = -v, veC\ 
c 

The Frechet derivative G^(0, 0) is an isomorphism on C^. Hence the Implicit Function Theorem 
applies and there exists a continuously differentiable map v : p ^ rj such that G{u{p)^p) = for all 
sufficiently small p. 

Next, we show that if p is the pressure consistent with the traveling water wave problem with 
velocity c, then ^{p) is indeed the corresponding water wave surface elevation. This is achieved 
by establishing that the reconstructed functions (j)^'^ and 4>^'^ are the horizontal and vertical fluid 
velocities (t)x and (j)z-, respectively. 

Let D = {— cxD < X < cxD, —h < z < r]}, as before, where rj represents the solution for the surface 
elevation of the traveling water wave problem with velocity potential (p. Hence <j) is harmonic in D. 
It is possible to harmonically extend (j) to D = {— oo < x < oo, —i] — 2/i < z < ry} by reflecting the 
problem across the mirror line z = —h. Thus (px is harmonic in D. 

If p is the pressure corresponding to the solution {(j), r]) of the water wave problem through 



(14), then at z = —h, (f)^'^ and its normal derivative take the same values as (j)x and its normal 
derivative, respectively. The Cauchy-Kowalevski Theorem for Laplace's equation [l4j implies that 
(px = (p^'^ ill a region near z = —h. But then (px and (f>^'^ and all their derivatives are equal in 
this region. This implies that (j)x = (p^'^ by analytic continuation. A similar argument shows 

that (pz = (p^'^ in D: to determine (pzz at z = —h we use the fact that due to the extension to D, 
(p is harmonic on z = —h. In addition, from (45) and (46), (p^'^ and (p^'^ harmonically extend up 



to z = M — h > ||r/||oo. Since (px = (p^'^ and (pz = (p^'^ in D, we can harmonically extend cpx and 
(pz Vl^ to z = M — h. This implies that (p^'^ = (px and cp^'^ = (pz at z = rj. Hence, from ([g]) in the 
traveling frame of reference. 



1 2 1 2 

-c(px + -(px + ^<Pz+9r] = 0, z = r] 

.ccpR,- + \{(P''n'' + 2 = 0, z = n, 



^G{r],p) = 0. 

From the Implicit Function Theorem, for small p all solutions r] to G{r],p) = are of the form 
rj = ^{p), where h' : p ^ rj is a map. Hence 



= (p'^'''{x,ri;p) = (px{x,r]) 



2gr] 



2gv{p) 



l + r]'x 



l + ul{p) 



where we used (13). In other words, there are functions rj which depend continuously on the true 
pressure p such that (21) is true. 

Next, assume there exists a different solution fj £ U C ,f] ^ r] such that 



4i¥ = 0«'(x,^;p), 
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for all X, where the pressure p is the pressure corresponding to the traveling water wave problem 
with velocity c. As before, (f>^'^ (x , fj; p) = (pxix^fj) and thus 



for all X. However, for a fixed c, the water wave problem has a unique traveling wave solution [3], 



which is contradicted by the statement that r] ^ fj. Thus the only solutions rj of (21) associated 
with the pressure p are the traveling wave solutions of the Euler equations. 



4 Asymptotic approximations 

In this section we derive a variety of asymptotic approximations to the pressure as a function of the 



surface elevation. Given the complexity of (21), such approximations are especially useful. In the 



sections below, we compare the results for the pressure obtained using (21), with those obtained 
from ([T]) and ([2]), as well as some asymptotic formulas obtained here. 
We introduce the nondimensional quantities z* , rj* and k*: 



C = C/L, z* = z/h, v* = r]/a, k* = Lk, c* = c/^gh, p*=p/gh, (47) 
where L is a typical horizontal length scale, and a is the amplitude of the surface wave. From (|3|-[6|), 



a nondimensional version of (21) is found to be 
/ — 2er] 1 



1 + (e/ir/5)2 2tt 



I e'^^Ti^Jc'^ - 2ep(e)| (A:) cosh (/uA; (1 + er/)) dk, (48) 

J — CO 



where e = a/h and /x = h/L. The *'s have been omitted to simplify the notation. This form of the 
nonlocal relation is our starting point to derive various approximate results. The two parameters e 
and jji provide many options for different asymptotic expansions: we may assume small amplitude 
waves (e <C //), or we may assume a long- wave approximation (// <C e), or we may balance both 
effects as in a Korteweg-de Vries (KdV)-type approximation (see [2], for instance). 

4.1 The small-amplitude approximation: SAOl and SA02 

If we expand r] in powers of e ^ 1, assuming that e ^ /i, we recover at leading order the approxi- 
mation 

r/(0 = cosh(Mfe)p(A)}} + O (e) , (49) 

where p{k) = J-{p}{k). Ignoring the O (e) term, ( |49[ ) is the nondimensional version of This 
demonstrates that ( |21[ ) is consistent with the frequently used ([2]), and we are able to recover such 
formulas in a consistent manner using the single equation (|21|) , instead of having to work with the 



full set of equations of motion. For the remainder of this paper, we will refer to the model (49) 
(solved for the pressure) as SAOl (Small- Amplitude, Order 1). 

If we proceed to higher order in e we find the presumably more accurate approximation 

v{i) = m{i) + ^m{i) + 0{e^), (so) 
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where 



m{0=J'''{cosHfik)m}}, (51) 

miO = - - + f^VoJ'-^kpik) smHfxk)} + ^T-^p^k) cosh(/iA:)}. (52) 



The formula (50) provides a new, exphcit, higher-order approximation for the surface elevation 
r/(^) in terms of the pressure p{(,) and the traveling wave speed c, assuming a small-amplitude 
approximation. For the remainder of this paper, we will refer to this model as SA02. 

4.2 The KdV Approximation: SWOl and SW03 

Alternatively, we can balance the parameters fi and e so that fi = yj~e. This is the KdV approxima- 
tion, see [HE])- At leading order, we recover the simplest approximation that the surface elevation 
equals the pressure: 

r/(e)=p(e) + 0(e). (53) 

This equation is exactly the hydrostatic approximation ([T]) in dimensionless variables; we will 
refer to this model as SWOl (Shallow Water, Order 1). Continuing the aproximations to higher 
order (up to order e^), we find 

''<«--ii-^(^i-i-KI)X=^-;^))-°<'''- 

We refer to this approximation as SW03 (Shallow Water, Order 3). 
Remarks. 



For G and p € Hm+s, it is possible to prove the analyticity of (48) in e and /x. Here, as 
before, M is related to the size of a symmetric strip around the real ^ axis. This analyticity 
serves to validate the asymptotic approximations derived above, as being obtained through a 
process that gives the first few terms of a convergent series. 



It appears to be a big restriction that the nonlocal formula (21) and the approximations 
derived above require a traveling wave profile. As shown in the next section, good results 
are also obtained for waves that are not merely traveling at constant speed. For waves in 
shallow water, excellent agreement is often obtained by using c = 1 (or c = \/g7i, returning 
to the dimensional version), which may be regarded as the zero-order approximation of an 
asymptotic series for c in terms of e. 

Using the procedures outlined in this section, the reader will find it straightforward to derive 
yet different approximations for the surface elevation in terms of the pressure measured at 
the bottom. For instance, one may consider a shallow- water approximation without imposing 
that the waves are of small amplitude, i.e., /x ^ e < 1, etc. 
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4.3 A heuristic formula: SA02h 



The transfer function approach ^ is very successful for a variety of reasons: (i) it is quite accurate, 
as is illustrated in the next few sections. This statement remains true to a varying degree for waves of 
relatively high amplitude; (ii) the most complicated aspect of using the formula is the computation 
of two Fourier transforms; and (iii) the formula applies to waves that are not necessarily traveling 
with constant speed. This is a consequence of the linearization that led to ([2]): each individual 
linear wave is traveling at constant speed, but typically their superposition is not. 

In this section we derive a different formula for the reconstruction of the surface elevation from 
the pressure at the bottom. This formula is obtained somewhat heuristically, and its justification 
rests on the fact that it agrees extremely well with both numerical and experimental data. Fur- 
thermore, its use requires the computation of only three Fourier transforms, and the velocity c does 
not appear in the final result. As a consequence, even though the derivation does not justify this, 
it is straightforward to apply to non-traveling wave profiles, where it performs very well. As for 
the other formulas above, the numerical and experimental results are presented below. 



An equivalent form of the nondimensional nonlocal equation ( 48 ) is 



1 - 2e7?/c2 1 



1 + (e/x7?g/c)2 2tt 



/oo 
e'''^P{k, e) cosh {fik (1 + erj)) dk, (55) 
-00 



where 



P{k,e) =t[i- Vl-2ep(e)/c2} (k). (56) 

So as to consider a small-amplitude approximation, we expand this equation in powers of e. How- 
ever, we do not expand P{k, e) at this point. Proceeding this way and retaining only first-order 
terms in erj and erj^, we find 



er/ 



er] 



2n 



e'-^^P{k,e) {cosh{fik) + e/ir//c sinh(^/c)) dk 



2^IZo^"''^Hk,e)cosh{fik)dk 



C2 271 



^ e*'=«P(A:, e)k smh{fik) dk 
Next, we expand P{k, e) in e, omitting terms of order and higher. We obtain 



(57) 



P{k,e) 



Substitution of (58) in (57) results in 



T] 



T {j>{k) cosh(;uA;)} 



(58) 



(59) 



1 — t^uT ^ {p{k) k smh{fik)} 

As stated above, this reconstruction formula does not depend on c, and its application requires 
the computation of a mere three Fourier transforms. This can be contrasted, for instance, with 
the formula SA02 which also uses a small-amplitude approximation. That formula requires the 
computation of five Fourier transforms and has explicit dependence on c. In fact, if one were to 



expand (59) in powers of e one would find at order the transfer function formula (|2|), and at 
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order the result SA02 with all c-dependent terms omitted. We refer to the results obtained 
using ((59]) as SA02h. 



5 Comparisons of the Different Approaches 

In this section, we present numerical results for the reconstruction of the surface elevation using 
the various relationships derived in Section [4] for both numerical and experimental pressure data. 
For the comparison using numerical data, we use previously computed periodic traveling waves 
solutions from [JJ^. By using the exact pressure underneath the traveling wave, we attempt to 
reconstruct the surface elevation. The same is done for various sets of experimental data obtained 
from the one- dimensional wave tank at the William Pritchard Fluids Laboratory at Penn State 
University. 



5.1 Comparison of the Different Approaches Using Numerical Data 

Using traveling wave solutions with periodic boundary conditions as calculated in [11], we determine 



the pressure at the bottom using (21) as follows. Without loss of generality, we assume that the 
solutions are 27r periodic. For a given traveling wave solution profile specified by {ijtrueiO ^ (^true) , 
the pressure p{(,) is obtained by equating the k-th Fourier coefficient of both the right- and left-hand 



side of (21) for k = —N . . . N, using a sufficiently high value of A^. This results in a linear system 
of algebraic equations for the coefficients of the Fourier series of \/ (? — 2p. Using this truncated 
Fourier series, we may solve directly for p{S,) in terms of the given solution set {rjtrueiO ■> ^^true) ■ Note 



that (25) offers a numerically equivalent alternative for computing p{C)- 

Our goal is to reconstruct the surface elevation from the thus computed pressure at the bottom, 
using the various formulas given above. The asymptotic formulas given in the previous section do 
not require anything more complicated than a fast Fourier transform. The solution of the nonlocal 



equation (21 ) is obtained using a pseudo-spectral method with differentiation carried out in Fourier 
space, while multiplication is carried out in physical space. We reconstruct by using a nonlinear 
solver such as a Gauss-Newton or Dogleg method [121 E2] with an error tolerance of 10~^^. As an 



initial guess for our nonlinear solver, we use the approximation from (50). Of course, the result 
obtained from the nonlocal equation should return the original surface elevation profile used to 
generate the pressure data, within machine precision. This provides a validation for the various 
numerical methods used. Here we compare the results from the asymptotic formulas of the previous 
section and evaluate their different errors. 

Using the parameter values /i = .l, g = l, /9 = 1 and L = 27r, we reconstruct the solution for 
various solution amplitudes and speeds. For solutions of small amplitude (say ak = .0001), we see 
that the reconstructions using all methods are in excellent agreement with the true surface wave 
elevation, see Figure [2^. However, even for waves with amplitudes less than 15% of the limiting 
wave height as given by [^ it becomes clear that certain approximations yield better results than 
others, see Figure ^p. In particular, while the nonlocal formula and the higher-order asymptotic 
formula reconstruct the wave profile well, the hydrostatic approximation SWOl reconstruction fails 
to reproduce an accurate reconstruction of the peak wave height. 

To demonstrate how the error changes as a function of the wave amplitude (or nonlinearity) , 
we compute the relative error 
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Figure 2: Reconstruction of the surface elevation from pressure data based on numerical experiments 
for h = 0.1, g = 1, p = 1 and L = 2tt. Amplitudes are a = 0.0001 (a) and a = 0.0056. (b). No 
legend is included: all approximations including the nonlocal formula (21 ) result in indistinguishable 
curves, except for the hydrostatic approximation SWOl, which displays a significant discrepancy 
for the bottom numerical experiment. 



error = — r| , (oOj 

I \ ^true \ |oo 

where Tjtrue represents the expected solution and r]r represents the reconstructed solution. For the 
same nondimensional parameters as before, we calculate the error as a function of increasing peak 
wave height demonstrated in Figure [3] As seen there, the error in all approximations grows as the 
amplitude of the Stokes wave increases. Figure |3] includes only solutions of small amplitude. If 
solutions of larger amplitude are considered, the discrepancies between the different approximations 
grow, as shown in Table |4| 

This table illustrates the large error generated by the lower-order methods SAOl and SWOl for 
waves which are no more than 55% of the limiting wave height as calculated in [7j . Even for waves 
which are 50% of the limiting wave height, the relative error (60) of the commonly used transfer 
function reconstruction SAOl exceeds 15%. In contrast, the higher-order methods SA02, SW03, 
and SA02h consistently yield more accurate results. It is also clear from the table that even for 
large amplitude waves, the nonlocal formula (21) (or more precisely for these numerical data sets, 
its periodic analogue (22)) provides a practical means to reconstruct the surface elevation from 
pressure data measured along the bottom of a fluid, at least in this numerical data setting. Below 
we establish the same using physical experiments. 

One limitation of the nonlocal equation (21) and some of its asymptotic counterparts from the 
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Calculated using c 
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Nonlocal Equation 

SW03 

SA01 

SA02h 




2.5 3 

Amplitude 



Figure 3: Plot of the relative error (60) in the reconstructed surface elevation rjj. as a function of 



the amplitude of rjtrue using the true value of the wave speed c. The asymptotic approximation 
SA02 is not included in this figure. It is more costly to compute than its heuristic counterpart 
SA02h, which yields better results. 



previous section is that they require the knowledge of the traveling wave speed c. In practice, this 
can be a difficult or impractical quantity to measure. One such impractical option is to include 
additional pressure sensors in order to measure the time it takes for the peak of the pressure data 
to travel from one sensor to another. A simpler option is to use approximations for the wave speed 
based on small-amplitude theory. For example, if we repeat the same error calculation as above, 
but with c ~ Vgh, we obtain Figure [s] We might hope that the reconstruction of the surface 
elevation would not suffer much. In fact, it appears unchanged. As seen in Figure [5} the error in 
reconstructing the peak wave height does not suffer at all from using this simple approximate (and 
amplitude-independent) value of c. Surprisingly, the error in the nonlocal reconstruction remains 
consistent with the error calculated using the true wave speed c. This lack of sensitivity to the 
precise value of c yields hope that with experimental data a simple approximation of the wave speed 
will be sufficient to accurately reconstruct the surface elevation. 



Percentage of 
Limiting Wave Height 


SWOl 


SAOl 


SW03 


SA02 


SA02h 


Nonlocal 


35 


21.14 


9.43 


4.01 


2.36 


1.76 


0.00 


45 


25.67 


13.43 


6.49 


4.27 


3.18 


0.00 


50 


27.88 


15.49 


7.89 


5.41 


4.04 


0.00 


53 


28.94 


16.49 


8.67 


5.96 


4.40 


0.00 


54 


29.65 


12.17 


9.15 


6.36 


4.70 


0.00 


55 


30.08 


17.58 


9.45 


6.61 


4.88 


0.00 



Table 4: Relative error (60) in percent, calculated comparing peak wave heights using various 



reconstruction formulas using numerical data. 
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Calculated using c. 




Figure 5: Plot of the error in the reconstructed surface elevation rjr as a function of the amplitude 
of rjtrue using an approximation for the wave-speed c. The SA02 model is not included for the 
same reason as in Figure [3j 



5.2 Comparison of the Different Approaches Using Data Prom Physical Exper- 
iments 



Here we discuss comparisons of results from the nonlocal formula (21) and the asymptotic ap- 
proaches with results from ten laboratory experiments performed at Penn State's Pritchard Fluid 
Mechanics Laboratory. In these experiments the pressure at the bottom of the fluid domain and 
the displacement of the air-water interface were measured simultaneously. The experimental fa- 
cility consisted of the wave channel and water, the wavemaker, bottom pressure transducers, and 
a surface displacement measurement system. The wavetank is 50 ft long, 10 in wide and 1ft deep. 
It is constructed of tempered glass. It was filled with tap water to a depth of h, as listed in 
Table [sj The pressure gage was a SENIORS PL6T submersible level transducer with a range of 
0-4 in. It provided a 0-5 V dc output, which was digitized with an NI PCI-6229 analog-to-digital 
converter using Lab View software. We calibrated this transducer by raising and lowering the water 
level in the channel. The pressure measurements had a high-frequency noise component, and thus 
were low-pass filtered at 20 Hz. The still-water height was measured with a Lory Type C point 
gage. The capacitance-type surface wave gage consisted of a coated-wire probe connected to an 
oscillator. The difference frequency between this oscillator and a fixed oscillator was read by a 
Field Programmable Gate Array (FPGA), NI PCI-7833R. Thus, no D/A conversion, filtering, or 
A/D conversion was required. The surface capacitance gage was held in a rack on wheels that are 
attached to a programmable belt. We calibrated the capacitance gage by traversing the rack at 
a known speed over a precisely machined, trapezoidally-shaped "speed bump". The waves were 
created with a horizontal, piston-like motion of a paddle made from a Teflon plate (0.5 in thick) 
inserted in the channel cross-section. The paddle was machined to fit the channel precisely with 
a thin lip around its periphery that served as a wiper with the channel's glass perimeter. This 
wiper prevented any measurable leakage around the paddle during an experiment. The paddle was 
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Figure 6: Wave tank comparisons of the reconstructed surface elevation with surface height mea- 
surements for h = 5.05 cm. This corresponds to Experiment #1 in Table [8] 

connected to the programmable belt and traveled in one direction. It was programmed with the 
horizontal velocity of a KdV soliton, which is given by 

u{x, t) = uq sech^ ( §—{x — Cot — uot/2)] , (61) 

V4/to Co J 

where cq = \/~gh, ao is the wave amplitude, and uq = aoco/h is the maximum horizontal velocity. 
The ao for the wavemaker displacement was varied between 2 cm and 3 cm. These values corre- 
sponded to large velocities and fluid displacements, outside of the regime of the KdV equation. The 
water adjusted to create a leading wave with a radiative tail. We compare results for the leading 
wave, where nonlinearity is likely to be important. 

To convert the time series of pressure and surface displacement data into spatial data, we 
use a combination of the sampling frequency and the estimated wave speed c. Specifically, let 
Pj represents the measured pressure at time tj = j At, where At is the time between pressure 
measurements. We assign a corresponding x value Xj to pj so that Xj = c{j At). From the 
pressure data measured from the physical experiments, we reconstruct the surface elevation using 
the same methods as in the previous section. For all experiments we use the admittedly simple 
approximation c ~ y/gh- We use the measured pressure data to reconstruct the surface elevation 



using the nonlocal formula (21), as well as the asymptotic approximations SWOl (hydrostatic), 
SAOl (transfer function), SA02, and SW03. 

As seen in Figure [6j the higher-order methods capture the peak wave height better than the 
lower-order methods, with the nonlocal equation yielding the most accurate representation of the 
peak wave height. The visual comparisons for all experiments is displayed in Figure[7} The nonlocal 



formula (21 ) consistently captures the peak wave height better than any of the approximate models 
derived in the previous section, and significantly better than the SWOl (hydrostatic) and SAOl 
(transfer function) models. This is quantified in Table ^ which displays the error (|60|) for the 



different approximations and the nonlocal equation (21). As is seen there, the result from (21) 



consistently produced the smallest error among all the reconstruction formulas. It is noteworthy 
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Figure 7: Wave tank comparisons of the reconstructed surface elevation with surface height mea- 
surements for various fluid depths. The experiments are ordered from left-to-right and top-to- 
bottom, corresponding to the experiment # in Table [8} 



that the heuristic SA02h approximation (59) consistently yields the second- lowest peak height 



error and consistently outperforms all other models except the nonlocal equation (21). Given the 



computational expense of solving the nonlocal equation, the approximation SA02h apparently 
yields the best compromise between efficiency and accuracy. 



6 Conclusion 



We have presented a new equation (21) relating the pressure at the bottom of the fluid to the 
surface elevation of a traveling wave solution of the one-dimensional Euler equations without ap- 
proximation. This equation is analyzed rigorously and the existence of solutions is proven using the 
Implicit Function Theorem. Solving the equation numerically is possible, but this is computation- 
ally relatively expensive when compared to currently-used approaches that require the computation 
of at most a few Fourier transforms. To this end, we derive various new approximate formulas, 
starting from the new nonlocal formula. The canonical approaches (hydrostatic approximation and 
transfer function approach) are easily obtained from the nonlocal formula as well. 



The different approximations and the nonlocal formula (21 ) are compared using numerical data, 
and their performance on physical laboratory data is examined. The nonlocal formula consistently 
outperforms its different approximations. For the numerical data this is by construction, as it was 
used to generate the numerical pressure data used for the comparison, starting from computed 
traveling wave solutions of the Euler equation. The higher-order approximate formulas result in a 
better reconstruction of the surface elevation compared to the hydrostatic or transfer function ap- 
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Experiment ^ 


Depth (cmj 


SWOl 


SAOl 


SW03 


SA02 


SA02h 


AT 1 1 

iNonlocai 


1 


5.05 


22.29 


6.51 


2.35 


0.84 


0.45 


0.20 


2 


5.05 


24.66 


8.05 


3.30 


1.13 


0.56 


0.01 


3 


5.05 


23.56 


7.75 


2.90 


0.95 


0.41 


0.18 


4 


5.05 


21.98 


7.18 


2.80 


1.29 


0.89 


0.66 


5 


5.05 


22.00 


6.63 


2.10 


0.48 


0.05 


0.03 


6 


3.55 


18.11 


5.21 


1.65 


0.65 


0.43 


0.36 


7 


3.55 


20.81 


7.04 


3.33 


2.22 


1.95 


1.52 


8 


4.10 


21.59 


8.39 


3.50 


2.18 


1.73 


0.93 


9 


4.10 


22.13 


7.65 


2.29 


0.33 


0.31 


0.05 


10 


4.10 


23.32 


9.25 


4.60 


2.91 


2.44 


2.13 



Table 8: Relative error (60) in percent, calculated comparing peak wave heights using various 



reconstruction formulas using experimental data. 



proaches. In the lab experiments, both the surface elevation and the bottom pressure are measured, 
allowing for an independent validation of the nonlocal equation (21). As expected, it outperforms 
the different approximations, where higher-order models perform better than lower-order ones. A 
good compromise between computational cost and obtained accuracy seems to be achieved by the 
heuristic approximation SA02h (59), which requires the computation of three Fourier transforms. 

Our derivation of the nonlocal equation ( 21 ) requires the surface elevation profile to be traveling 
with constant speed c. Regardless, we show that the results are not sensitive to the exact value of 
c and even rough estimates (i.e., c = \/gK) provide excellent results, both for the nonlocal equation 



(21) and its various asymptotic approximations, most notably (59). 
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